Probing white dwarf interiors with LISA: periastron precession in double white dwarfs 
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In globular clusters, dynamical interactions give rise to a population of eccentric double white 
dwarfs detectable by the Laser Interferometer Space Antenna (LISA) up to the Large Magellanic 
Cloud. In this Letter, we explore the detectability of periastron precession in these systems with 
LISA. Unlike previous investigations, we consider contributions due to tidal and rotational dis- 
tortions of the binary components in addition to general relativistic contributions to the periastron 
precession. At orbital frequencies above a few mHz, we find that tides and stellar rotation dominate, 
opening up a possibly unique window to the study of the interior and structure of white dwarfs. 

PACS numbers: 95.10.Ce, 95.30.Sf, 95.55.Ym, 97.10.Cv, 97.60.-s, 97.80.-d, 



Introduction - Binary stars have long been recognized 
as unique astrophysical laboratories for the study of 
physics and cosmology [l| . When one of the binary com- 
ponents is a compact object, measurements of its inter- 
actions with the orbital companion provide a wealth of 
information on the compact object and the properties of 
matter under extreme conditions. 

The Laser Interferometer Space antenna (LISA; 0) 
will survey the whole galactic population of binaries con- 
sisting of two white dwarfs (WDs) with periods ^ 6 hi 
by monitoring their gravitational wave (GW) emission. 
LISA will individually resolve ^ double white dwarfs 
(DWDs) 3] and will be particularly effective in discov- 
ering short-period ( ^ 30 min) systems - in stark con- 
trast to current and planned electromagnetic observato- 
ries. LISA will therefore allow astrophysical studies of 
outstanding questions in compact object binary physics, 
such as dynamically unstable mass transfer, accretion, 
and type la supernova formation [^] . 

So far, investigations of DWDs have focused on the 
formation and properties of binaries with circular orbits, 
which dominate the Galactic DWD population. How- 
ever, Willems ct al. jH recently predicted the existence 
of a sub-population of DWDs consisting of eccentric bi- 
naries formed through dynamical interactions in globu- 
lar clusters. Unlike present and planned electromagnetic 
telescopes, LISA will be able to detect these systems to 
distances as far as the Large Magellanic Cloud, providing 
a unique opportunity to study degenerate matter through 
the imprint of tidal effects in the detected GW signal. 

Tidal interactions in close binaries couple the spins of 
the component stars to the orbital motion, driving bi- 
naries to a state of minimum kinetic energy in which 
the orbit is circular, the stellar spin angular momenta 
are aligned with the orbital angular momentum, and the 
stellar rotation rates are synchronized with the orbital 
motion. The efficiency of this process depends strongly 
on the mode of energy dissipation in the stellar interiors, 
which is reasonably well understood for non-degenerate 
stars, but an open question for WDs However, in ec- 
centric binaries, tides also cause a non-dissipative preces- 
sion of the periastron of the orbit. This "apsidal motion" 
is caused by perturbations in the gravitational field due 



to the tidal distortions of the binary components. The 
apsidal motion depends primarily on the internal mass 
distribution of the stars, and takes place on considerably 
shorter time scales than dissipative spin-orbit coupling. 
The shorter time scales facilitate the inference of apsidal- 
motion rates from electromagnetic and GW observations. 

In addition to tides, rotational distortions of binary 
stars and general relativity (GR) also contribute to the 
apsidal motion in eccentric binaries. While the rotational 
contribution depends on the internal mass distribution 
in a similar way as the tidal effect, the GR contribution 
depends only on the total system mass and orbital el- 
ements. For non-degenerate stars, comparisons between 
theoretically predicted apsidal-motion rates and observa- 
tionally inferred rates have long served as a critical test 
of theories of stellar structure and evolution 0, S ■ 

In this Letter, we examine the physics accessible by 
measuring the apsidal motion of eccentric DWDs with 
LISA. General relativistic apsidal motion has already 
been proposed as a tool to derive the total system mass 
of eccentric neutron star (NS) binaries [ll|. Applications 
to WD binaries have so far not been considered, nor have 
the imprint on GWs of the periapse precession induced 
by tides and stellar rotation. 

Apsidal motion. — We briefly summarize the equations 
governing the tidal, rotational, and GR contributions to 
the apsidal motion in eccentric binaries. For this purpose, 
we consider a binary consisting of two uniformly rotating 
stars with masses Mi. 2, radii -R1.2, and rotational angular 
velocities J7i.2.The stellar rotation axes are assumed to be 
perpendicular to the orbital plane [13] ■ We furthermore 
let P be the orbital period, a the semi-major axis, e the 
orbital eccentricity, and 7 the argument of the periastron. 

The contribution to the apsidal motion from the tidal 
distortion of the binary components is most commonly 
determined under the assumption that the orbital and 
rotational periods are long compared to the periods of the 
free oscillation modes of the component stars [l2„ J3,, jj] . 
Under this assumption, the rate of secular apsidal motion 
due to the dominant quadrupole tides raised in star i 
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(i = 1, 2) is given by 
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where ki is the quadrupolar apsidal-motion constant of 
star i. When the orbital and rotational periods are of the 
order of the free oscillation modes of the binary compo- 
nents, deviations from Eq. ([T]) arise due to the increasing 
role of stellar compressibility on the tidal displacement 
field for higher tidal forcing frequencies and due to the 
occurrance of resonances between dynamic tides and non- 
radial stellar oscillations [l^ . 

The apsidal-motion constants ki measure the degree to 
which mass is concentrated towards the stellar center and 
are determined by numerical integration of the equation 
of Clairaut (for details see, e.g., [3, HI])- The constants 
are unaffected by dissipative effects as long as the tidal 
forcing is not in resonance with any of the WDs' non- 
radial stellar oscillation modes 1^, 21 1. In the limiting 
cases where the stars are approximated by point masses 
or equilibrium spheres with uniform mass density, the 
constants ki take the values and 0.75, respectively. For 
more realistic stellar models, the constants take values 
between these two extremes. 

Rotation contributes to the apsidal motion through the 
rotational quadrupole distortion caused by the centrifu- 
gal force [ij]. The corresponding rate of secular apsidal 
motion depends on Ri^2 and fci_2 in a similar way as the 
tidal contribution to the apsidal-motion rate, but has a 
different dependence on Mi. 2 and e: 
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where = 2 tt/P is the mean motion. 

The GR contribution to the apsidal-motion rate differs 
from the tidal and rotational contributions in that it is 
independent of the radii and internal structure of the 
binary components. At the leading quadrupole order, 
the GR apsidal-motion rate is given by 



7gr = 



2 TT 3 G + M2 
"P ~ a(l-e2)' 
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where G is the Newtonian gravitational constant, and c 
the speed of light [T^]. 

In Fig. [H the tidal, rotational, GR, and total apsidal- 
motion rate are shown as functions of the orbital fre- 
quency V = \/P for different orbital eccentricities and 
conservative WD component masses of 0.3 and 0.6 Mq. 
The tidal and rotational contributions are determined 
using Nauenberg's j20| zero-temperature mass-radius re- 
lation, and setting /ci = A:2 = 0.1. These fci and ^2 are 
appropriate for cool WD models of 0.3 and 0.6 Mq; 
the dependence of ki and ^2 on the WD mass and tem- 
perature will be explored in more detail in a separate 
investigation. The WD rotational angular velocities are 
furthermore assumed to be synchronized with the orbital 
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FIG. 1: Left: Tidal (short-dashed), rotational (dotted), GR 
(long-dashed), and total (solid) apsidal-motion rate for DWDs 
with Ml = M2 = 0.3 M0. Right: Total apsidal-motion rate 
for DWDs with Mi = M2 = 0.3 Mq (solid) and Mi = M2 = 
0.6 Af© (short-dashed). Both panels show results for e = 0.01 
(black) and e = 0.5 (grey), and ki = k2 = 0.1. 



angular velocity at periastron. Since the tidal effects usu- 
ally dominate the rotational effects, this assumption does 
not affect the main conclusions of the calculation. 

It is evident that the total apsidal-motion rate is sub- 
stantial throughout the entire LISA band, even for eccen- 
tricities as low as e~ 0.01: ati^~0.5 mHz, periastron 
precession already induces a phase shift in the GW sig- 
nal of more than 27r over an observation time Tobs = 5 yr 
(the current minimum mission lifetime requirement) and 
therefore becomes potentially detectable. The phase shift 
is even larger at higher orbital frequencies. Equally strik- 
ing is the dominance of tides and stellar rotation at fre- 
quencies above a few mHz. Willems et al. [Bj have shown 
eccentric DWDs in this frequency range to be detectable 
by LISA up to distances as far as the Large Magellanic 
Cloud, opening up new avenues for GW astrophysics of 
DWDs. The tidal and rotational contributions to 7 fur- 
thermore decrease with increasing mass of the WDs due 
to the smaller radii of more massive WDs. At low fre- 
quencies, where GR effects dominate, the apsidal-motion 
rate increases with increasing mass of the WDs. 

LISA observations. — Periastron precession leaves a 
signature in the GW forms of eccentric DWDs by mod- 
ifying the phase of the signal recorded by laser interfer- 
ometers. LISA can therefore probe into the structure of 
DWDs by observing the apsidal motion due to the above 
contributions. Here, we first show conceptually how one 
can measure 7 and the parameters that drive its evolu- 
tion; next we explore how accurately these parameters 
can be measured. Throughout this discussion we model 
GW radiation at the leading Newtonian quadrupole or- 
der (post-Newtonian corrections are negligible in this fre- 
quency and mass range), and model LISA following p^ . 

The signal from an eccentric DWD in the LISA frame 
can be schematically written as 
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where J+,x sjce the antenna beam patterns (that depend 
on the source right ascension a and decUnation d, and 
the wave polarization ip sd, a, constant reference time), 
and 

= A|-(l + cos^i)ii„(e)cos[ri.(/.(i) + 27(t)] 
— (1 + cos^ t)wri(e) cos[n(/)(i) — 2'-f{t)] 
+ sin^ tw„(e) cos[n0(t)]|, (5) 

h':^^{t) = 2Acosi^Un{e)sm[n(j){t) + 2-/{t)] 

+w„(e)sin([n(/)(t)-27(i)])}. (6) 



(f){t) is the Dopplcr shifted orbital phase 
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Here, 

2TTi/t + TTiyt^ + O {t^) + (po , where i' is the orbital frequency 
derivative and (j)o an arbitrary initial phase; t is the con- 
stant source inclination angle, A = (2'!:v)^l'^ j d the 
GW amplitude, the chirp mass, and d the distance 
to the source; u„(e), Un(e), and i«„(e) arc linear combi- 
nations of the Bessel functions of the first kind J„(ne), 
Jn±iine) and Jn±2(^^e); explicit expressions can be de- 
rived using Eqs. (7) and (10) in [l6l |. 

In the absence of periastron precession, radiation is 
emitted at multiples n of the orbital frequency but 
when periastron precession is present, each of these har- 
monics is split into a triplet with frequencies niy±'j/TT and 
nv, and amplitudes u„(e), u„(e) and ^^(e), respectively. 
As already noted by the observation of any two 

"emission lines" , allows us to derive the orbital and the 
apsidal-motion frequency. In practice, for typical DWD 
eccentricities (e ^ 0.5; see @), |'«n(e)| ^ |wn(e)| , |?i'„(e)|, 
so that LISA will primarily rely on observations of GWs 
at frequencies + j/ir for at least two values of n. 

We can compute whether and how accurately perias- 
tron advance can be measured by computing the Fisher 
information matrix (see e.g. jig ) associated with the 
measurement. The signal depends on a, S, l, ip, A, e, 
V, i/, 7, (j)Q and 70 (the argument of periastron at an ar- 
bitrary reference time). Due to the fact that a, S, l and 
■0 are only weakly correlated with the remaining param- 
eters for an observation lasting several years |18| . we do 
not include them in our analysis and compute the angle- 
averaged Fisher information matrix for Tobs = 5 yr. We 
conservatively consider the first 10 harmonics in Eq. ^ 
and normalize the results to signals detected at an opti- 
mal signal-to-noise ratio (SNR) of 10 (the mean-square 
errors scale as « 1/SNR). The results are summarized in 
Table U for the parameters relevant to this investigation. 

For v ^ 1 mHz, the apsidal motion becomes progres- 
sively harder to measure due to the smaller and smaller 
phase shift (cf. Fig.[T]). At w 0.1 mHz the estimated er- 
ror on -j/n is greater than 7/77 itself and apsidal motion 
becomes undetectable (the details of course depend on 
the actual values of e, z^. Mi, M2 and SNR for the source 
at hand) [2^. However, for 1/ ^ ImHz, LISA can detect 
periastron precession and measure 7/7r with a relative 
error of ~ 1% — 10%, depending on the mass and ec- 
centricity of the binary. This result holds even for small 
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TABLE I: Measurement accuracy of e, "//n, and u for selected 
values of e and v. The errors are for Tobs = 5yr and for a 
source detected with optimal SNR = 10. No assumption is 
made on the physical mechanism driving the apsidal motion 
7 and/or frequency drift v. The error on v is reported only 
for systems with v > 1 mHz, where the GR contribution to !> 
becomes observable. 
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TABLE H; Relative error on total system mass and chirp mass 
measurements assuming only GR contributes to 7 and !>, for 
a 5-yr observation of a Mi^M2 = 0.6Mq DWD at SNR = 10. 



eccentricities e ~ 0.01 At these frequencies most 

systems are also expected to exhibit a detectable change 
of the orbital frequency. Assuming that only general rel- 
ativity affects 7 and i/, the combined measurement of 7 
and v allows us to determine the total system mass M 
and chirp mass M , as shown in Table |TT1 and therefore 
the individual WD masses Mi and M2. 

Astrophysical implications. — We have investigated the 
impact of tides, rotation, and GR on the apsidal motion 
of eccentric DWDs and its signature on the emitted GWs. 
Based on our present understanding of the astrophysical 
scenarios @, we conclude that LISA will be able to ob- 
serve the periastron advance for the vast majority of such 
sources detected. These observations provide a new and 
unique probe into the internal structure of WDs. 

Tides and stellar rotation strongly dominate the 
apsidal-motion rate at orbital frequencies above ~ 1 
mHz, inducing phase shifts much larger than those esti- 
mated using only the GR contribution. In GW searches 
for eccentric binaries, it is therefore essential to include 
7 in the signal templates as a phenomenological param- 
eter not bound by 7gr in order to not bias LISA surveys 
against eccentric DWDs. In the interpretation of the 
data, neglecting tidal and rotational contributions would 
lead to an overestimate of the total system mass derived 
from 7. This will likely induce a misclassification of GW 
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sources as NS rather than DWD binaries, thus affecting 
the ratios of different populations of compact object bi- 
naries that hold essential signatures of stellar evolution 
and binary formation mechanisms. On the other hand, 
the dependence of 7 on i?i,2, -^^1,2, and fci^2 provides a 
new window into the internal structure of WDs and the 
astrophysics of such stars, but poses a severe degeneracy 
problem for the extraction of astrophysical information 
from measured apsidal-motion rates. More refined theo- 
retical modeling is therefore needed in preparation of the 
LISA mission to fully characterize the dependence of 7 
on the WD physical parameters and to identify routes 
to untangle them. Even though we have here focused on 
DWDs, our results also apply to WDs with NS compan- 
ions. These sources are in fact much "cleaner" probes of 
WD physics since the tidal and rotational distortions of 
the NS contribute negligibly to the apsidal motion. In 
this case, the apsidal-motion rate therefore carries the 
unique signature of only one WD rather than two. 

At frequencies ^ 1 mHz, GR effects dominate the pre- 
cession rate and measurements of 7 allow the determi- 



nation of the total system mass [TT|. For 1/ ^ 0.5 mHz 
0, the radiation reaction may also cause a measurable 
drift in the orbital frequency. Assuming that there is 
no significant contribution from tidal and/or magnetic 
spin-orbit coupling, measuring 1) with LISA yields the 
source distance as well as the chirp mass [e.g. ha. In the 
(small) regime where general relativity dominates the ap- 
sidal motion and ij is detectable, the combined knowledge 
of the chirp mass and total system mass yields the masses 
of the individual WD components. For eccentric NS-NS 
binaries - which have negligible tidal/rotational distor- 
tions and tidal/magnetic spin-orbit coupling throughout 
the LISA band and a stronger z> than WD binaries - the 
measurement of the individual masses at w 10% level 
should be routine. 
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